Anderson localization of a Tonks-Girardeau gas in potentials with controlled disorder 



O 
(N 



(N 



C/2 



I 

o 



^, 

o 
o 



X 



J. Radic/ V. Bacic/ D. Jukic/ M. Scgev,^ and H. Buljan^ 

'Department of Physics, University of Zagreb, PP 332, 10000 Zagreb, Croatia 

^Techmon, Israel Institute of Technology, Haifa, Israel 

(Dated: May 27, 2010) 

We theoretically demonstrate features of Anderson localization in the Tonks-Girardeau gas con- 
fined in one-dimensional (ID) potentials with controlled disorder. That is, we investigate the evo- 
lution of the single particle density and correlations of a Tonks-Girardeau wave packet in such 
disordered potentials. The wave packet is initially trapped, the trap is suddenly turned off, and 
after some time the system evolves into a localized steady state due to Anderson localization. The 
density tails of the steady state decay exponentially, while the coherence in these tails increases. 
The latter phenomenon corresponds to the same effect found in incoherent optical solitons. 

PACS numbers: 05.30.-d, 03.75.Kk, 67.85.De 



I. INTRODUCTION 

The phenomenon of Anderson localization |l|, which 
was originally theoretically predicted in the context 
of condensed matter physics, has been experimentally 
demonstrated in other wave systems including optical 
waves [2ny] and ultracold atomic gases (matter waves) 
[3, Q. In the context of Bose- Einstein condensates 
(BECs), Anderson localization was obtained by plac- 
ing ultracold atomic BECs in elongated, essentially one- 
dimensional disordered (7[ and quasiperiodic incommen- 
surate potentials Q, which were created optically (see 
Ref. Q for a recent review of the topic). The mat- 
ter waves utilized in those experiments were conden- 
sates, i.e., they were spatially coherent in the sense 
that their one-body density matrix factorizes p{xi, X2) ~ 
$*(xi)$(a;2), where $(x) is the condensate wave func- 
tion. However, in reality interactions and/or the presence 
of the thermal cloud affects the spatial coherence in the 
system. Naturally, the spatial coherence in the system 
is expected to have important implications on localiza- 
tion phenomena, since the phenomenon of Anderson lo- 
calization is deeply connected to interference of multiple 
reflected waves. This motivates us to study Anderson lo- 
calization in a Tonks-Girardeau gas, which is a relatively 
simple example of partially-spatially-coherent Bose gas 
(i.e., it is not condensed). 

The Tonks-Girardeau model describes a system of 
strongly repulsive (" impenetrable" ) bosons, confined in 
one-dimensional (ID) geometry [l3|- Exact solutions of 
the model are found by employing the Fermi-Bose map- 
ping [id . Ill| . wherein the Tonks-Girardeau wave func- 
tion (for both the stationary and the time-dependent 
problems) is constructed from a wave function describ- 
ing noninteracting spinless fermions. In Ref. |12| it was 
suggested that the Tonks-Girardeau model can be exper- 
imentally realized with ultracold atoms in effectively ID 
atomic waveguides. This regime is reached at low tem- 
peratures, for sufficiently tight transverse confinement, 
and with strong effective interactions [13 - [l4| . Indeed, in 
2004 two groups have experimentally realized the Tonks- 



Girardeau gas |15l.ll6l|. Furthermore, nonequilibrium dy- 
namics of a ID Bose gas (including the Tonks-Girardeau 
regime) has been experimen tally addressed in the context 
of relaxation to equilibrium [17J | . It is known that ground 
states of the Tonks-Girardeau gas on the ring [1^, or in a 
harmonic potential [l^ are not condensates, because the 
population of the leading natural orbital scales as VN, 
where N is the number of particles. Thus, the Tonks- 
Girardeau gas is only partially spatially coherent. The 
free expansion of the Tonks-Girardeau gas from some ini- 
tial state has been of great interest over the past few years 
[20I - I23I I : this type of scenario, i.e., expansion from an ini- 
tial state which is localized (say by a trapping potential) 
can be used to address Anderson localization [7[. 

The experimental demonstrations of Anderson local- 
ization in ultracold atomic gases were preceded bytheo- 
retical investigations of this topic (e.g., see Refs. |24l - [26| . 
see also Ref. [9| and references therein) . The interplay of 
disorder (or quasiperiodicity) and interactions in a Bose 
gas (from weakly up to strongly correlated regimes) , has 
been often studied in the context of the Bose-Hubbard 
model [H m, [27l - [37t . Within the model, a transition 
from a superfluid to a Bose glass phase has been pre- 
dicted to occur [23, [2^1 . The aforementioned interplay 
has been studied by using versatile methods including 
calculating the energy absorption rate [37[, momentum 
distribution and correlations [30, ISa] , and expansion dy- 
namics [33|, yj]. In the limit of strong repulsion, the 
system can be described by using hard-core bosons on 
the lattice [30, S, [33] ■ For these systems, by employing 
the Jordan- Wigner transformation the bosonic system is 
mapped to that of noninteracting spinless fermions, and 
all one-body observables can be furnished from the one 
body density matrix both in the stationary (e.g., see [SSj ) 
and out-of-equilibrium systems [20|. The ground state 
properties of the hard-core Bose gas in a random lattice 
have been studied in [SOf, whereas expansion dynamics 
was considered in |33|; both approaches predict the loss 
of quasi long-range order. 

Here we study Anderson localization within the 
framework of the Tonks-Girardeau model \X^ in one- 
dimensional disordered potentials. We study the expan- 



sion of a Tonks- Girardeau wave packet in a potential with 
controlled disorder. The potential is characterized by its 
correlation distance parameter a. Pd t — 0, the initial 
wave packet is in the ground state of a harmonic trap with 
frequency w (with small disorder superimposed upon it), 
and then the trap is suddenly turned off. After some 
time, we find that the system reaches a steady state char- 
acterized by exponentially decaying tails of the density. 
We show that the exponents decrease with the increase 
of w and the decrease of a in the investigated parameter 
span (cr = 0.13 - 0.40 /zm and w = 5 - 10 Hz). The 
one-body density matrix pB{x,y,t) of the steady state, 
that is its amplitude \pB{0,x,t)\, decays exponentially 
on the tails of the localized wave packet. However, in 
the region of these tails the degree of first order coher- 
ence |/iB(0, a;,i)| = \pb{0, x, t)|/^/pB(0, 0, t)pB{x, x, t) 
reaches a plateau. These plateaus are connected to the 
behavior of the single-particle states used to construct the 
Tonks-Girardeau wave function, from which we find that 
the spatial coherence increases in the tails. This increase 
of coherence in the tails has its counterpart in incoherent 
optical solitons [39[, a phenomenon well understood in 
terms of the modal theory for incoherent light [S^ • 

The paper is organized as follows: In Section |ll] we 
describe the Tonks-Girardeau model. In Section IIIII we 
present our numerical results on Anderson localization 
in this system. Finally, in Section IIVI we outline our 
conclusions. 



The solution of this system may be written in compact 
form via the famous Fermi-Bose mapping, which relates 
the Tonks-Girardeau bosonic wave function ij^b to an an- 
tisymmetric many-body wave function ij^p describing a 
system of noninteracting spinless fermions in ID [lO[: 



-0B = Iil<i<j<NSgn{xi - Xj)iPf{xi,X2, ■ ■ .,XN,t). (2) 

In many physically relevant situations, the fermionic 
wave function ipp can be written in a form of the Slater 
determinant. 
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where ipmix^t) denote N orthonormal single-particle 
wave functions obeying a set of uncoupled single-particle 
Schrodinger equations 
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i^^{x,t), m = l,...,iV. (4) 



In Eqs. ©-(HI) we have outlined the construction of 
the many-body wave function describing the Tonks- 
Girardeau gas in an external potential V{x)^ both in the 
static [lO] and time-dependent case (llj . 

Given the wave function ipB, we can straightforwardly 
calculate all one-body observables furnished by the re- 
duced single-particle density matrix (RSPDM), 



II. TONKS-GIRARDEAU MODEL 

In this section we present the Tonks-Girardeau model 
which describes "impenetrable-core" ID Bose gas [10, 
[lll |. We study a system of A^ identical Bose particles 
in ID geometry, which experience an external potential 
V{x). The bosons interact with impenetrable pointlike 
interactions [lO[ , which means that the wave function de- 
scribing the bosons vanishes whenever two particles are 
in contact, that is, 4'b{xitX2t ■ ■ -.XNTt) = if Xi = Xj 
for any i ^ j. The wave function ij^b must also obey the 
Schrodinger equation 
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here wc use dimensionless units as in Ref. |40| . i.e., x = 
X/Xo, t = T/Ta, and V{x) = U{X)/Eo, where X and T 
arc space and time variables in physical units, and U{X) 
is the potential in physical units. Given the particle 
mass 771, the time scale Tq = 2mXQ/h and energy-scale 
Eq = fi^ / {2mXQ) are set by choosing an arbitrary spatial 
length-scale Xq. For example, in our calculations below 
Xo = 1 pm, while the mass corresponds to *^Rb, which 
yields the temporal scale Tq = 2.76 x lO"'^ s, and the en- 
ergy scale £^0 = 3.82 x 10""^^ J. The wave functions are 
normalized as J dxi . . . dxN \ ipB {xi , a;2 , . . . , xn , i) ^ = 1 . 



PBix,y,t) 



N dx2 . ■ . dxN i'six, X2,---, xn, t)* 

X1pB{y,X2,---,XN,t); (5) 



by employing the formalism presented in Ref. |4l| . If the 
RSPDM is expressed in terms of the single-particle wave 
functions ipm as 



N 



PB{x,y,t)^ ^ ip*{x,t)A^j{x,y,t)il;j{y,t), (6) 



it can be shown that the N x N matrix A{x,y,t) 
{Aij{x,y,t)} has the form 



A{x,y,t) = (P"^)-* detP 



(7) 



where the entries of the matrix P are Pij{x,y,t) = 
dij — 2 [^ dx'^*(x\ t)ipj{x', t) (x < y without loss of gen- 
erality) [U . 



III. NUMERICAL RESULTS ON ANDERSON 
LOCALIZATION IN A TONKS-GIRARDEAU 

GAS 

In order to investigate Anderson localization of the 
Tonks-Girardeau gas, we perform numerical simulations 
designed in the fashion of optical [5[ and matter wave Q 



experiments which were conducted recently to demon- 
strate Anderson localization. We investigate dynamics 
of a Tonks-Girardeau wave packet in a disordered poten- 
tial Voi'x), where the initial wave packet (at i = 0) is 
localized in space by some trapping potential. After long 
time of propagation the wave packet reaches some steady 
state. Anderson localization is indicated by the exponen- 
tial decay of the density of the wave packet in this steady 
state. 

More specifically, we assume that initially, at i = 0, 
the gas is in the ground state of the harmonic oscillator 
potential, with the small controlled disordered potential 
superimposed upon it, that is, 



V{x) = Vd{x) + v^x^ for t < 



(8) 



At t = the trapping potential is suddenly turned off, 
i.e., 



V{x) = Vd{x) for t > 0, 



(9) 



after which the density and correlations of the gas begin 
to evolve. This means that at i = the wave function ips 
is given by Eqs. ^ and ^ where ^pm{x,t = 0) is the 
mth single-particle eigenstate of the potential Vd{x) + 
v^x^ . The subsequent evolution of V'm is given by Eq. 
dl]) where the potential is given solely by the disordered 
term V{x) = Vd{x). 

The construction of the disordered potential Vd{x) 
used in our numerical simulations is described in the Ap- 
pendix. The disordered potential can be characterized 
in terms of its correlation functions; the autocorrelation 
function is defined by 



Acix) = {VD{x'-x)VDix')), 



(10) 



where Vd{x) = Vd{x) — (^0(2;'))^'' ^^^^l (• • • )xi denotes 
a spatial average over x' . For the disordered potentials 
in our simulations we have approximately 



Ac{x) = Vo' 



2 sin {x/ a) 



{x/a) 



(11) 



where a denotes the spatial correlation length of the dis- 

2 

ordered potential, whereas Vq = {V ij(x))x denotes its 
amplitude. The spatial power spectrum of the potential 
has support in the interval [—Kcut^ ^c«t]j where the cut- 
off value is Kcut = 'i^/cf. Thus, the potential Vd{x) has 
the autocorrelation function identical to that of the op- 
tical speckle potentials used in the experiments, e.g., see 
0. 

The asymptotic steady state of the system depends on 
the parameters of the disordered potential a and Vq , and 
on the initial state, that is, the harmonic trap parameter 
v. In fact, since the dynamics of the Tonks-Girardeau 
gas is governed by a set of uncoupled Schrodinger equa- 
tions, it follows from the simple scaling of units outlined 
below Eq. ([T]), that there are in fact only two indepen- 
dent parameters; thus we investigate the dynamics in de- 
pendence of v and cr, and keep Vq at an approximately 
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FIG. 1: (color online) Anderson localization in a Tonks- 
Girardeau gas in dependence of the initial trap parameter v 
(i.e., uj). The parameters of the disordered potential are (a = 
0.13 and Vb = 0.465). (a) The averaged density of the Tonks- 
Girardeau wavepacket at f = and after t — 1450 (=4 s) of 
propagation. The initial state corresponds to :/ = 8.67 x 10~^ 
(id — 10 Hz), (b) Shown is the density of a Tonks-Girardeau 
gas (in the localized steady state) after t = 1450 (=4 s) of 
propagation in a disordered potential. Blue dot-dashed line 
corresponds to ;/ = 4.34 x 10"'^ {lj = 5 Hz), whereas red solid 
line corresponds to :/ = 8.67 x 10~^ {uj = 10 Hz), (c) Same 
as figure (b) on a logarithmic scale. Blue line corresponds to 
{uj = 5 Hz), and red line corresponds to u = 10 Hz; arrow 
indicates the increase of oj. For |a;| larger than some value 
(call it Lt), the density decays exponentially, which charac- 
terizes Anderson localization. The density-tails decay slower 
for larger initial trap parameter u; (see text for details). 



constant value. We have performed our numerical sim- 
ulations in a region of the parameter space which was 
accessible with our numerical capabilities, but which is 
relevant to experiments [3, 13 . We have varied the cor- 
relation length a of the potential from 0.13 up to 0.40 




FIG. 2: (color online) Anderson localization in a Tonks- 
Girardeau gas in dependence of the disorder parameter a. 
The averaged density of a Tonks-Girardeau gas after t — 1450 
(=4 s) of propagation in the disordered potential. The plots 
correspond to cr = 0.13 (blue line), a = 0.25 (green line), and 
a = 0.40 (red line); arrow indicates the increase of a. The 
initial state corresponds to a; = 10 Hz, while the amplitude 
of the disordered potential is Vb ~ 0.47 (see text for details). 
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FIG. 3: (color online) The exponential decay rates 7'^^ (dot- 
ted blue line) and 72^+7 (solid red line), obtained with per- 
turbation theory [42l, in dependence of the correlation length 
of the potential a. Circles represent decay rates obtained nu- 
merically. The initial trap parameter is cj = 10 Hz (see text 
for details). 



(corresponding to 0.13 /im and 0.40 /im since the spatial 
scale is chosen to be Xq = 1 /im) , and the harmonic trap 
parameters in the interval ly = 4.34 — 8.67 x 10~^ (cor- 
responding to cj = 5 — 10 Hz). The number of particles 
used in our simulations is relatively small, N — 13, due 
to the computer limitations, however, despite of this, one 
can use our simulations to infer general conclusions that 
would be valid in an experiment with larger N. It should 
be emphasized that all plots of densities and correlations 
are ensemble averages made over 40 realizations of the 
disordered potentials. 

First we investigate the behavior of the single-particle 
density. In Figure [1] we show pB{x,x,t = 1450) ver- 
sus X for (cr, Vo) = (0.13,0.465), and two values of v. 
V ^ 4.34 X 10^2 (w = 5 Hz) and i^ = 8.67 x 10"^ 
{u! = 10 Hz). In Fig. [TJa) wc compare the initial density 
(at t = 0), with the density at time t = 1450 (= 4 s). 



at which the steady state regime is already achieved (all 
graphs below which describe the steady state are also cal- 
culated at this time). We observe that the steady-state 
density has a broad central part with a fairly flat top, and 
decaying tails on its sides. The central part is composed 
of many single-particle states ipj . In Fig. [TJb) we plot the 
steady state density for two values of uj. For oj = 5 Hz, 
the central part is broader than for a; = 10 Hz, but the 
tails are decaying faster with the increase of |a;|, as shown 
in Fig. IHc), where the densities are plotted in the loga- 
rithmic scale. We clearly see that for |a;| larger than some 
value (call it Lj), the density decays exponentially, which 
indicates Anderson localization. We have fitted the tails 
to the exponential curve pB{x,x,t) ex exp(— A|a;|) and 
obtained A = 0.0097 for w = 5 Hz, and A = 0.0053 for 
a; = 10 Hz, that is, we find that the density-tails de- 
cay slower for larger initial trap parameter lu. For larger 
values of uj, the trap is tighter and the initial state has 
larger energy and broader momentum distribution, there- 
fore, it is harder to achieve localization of the wave packet 
(e.g., see [2^|42]). Another way to interpret these sim- 
ulations is in terms of the spatial correlation distance of 
the wave packet. An incoherent wave packet can be char- 
acterized by using the spatial correlation distance, which 
determines a spatial degree of coherence; this quantity is 
inversely proportional to the width of the spatial power 
spectrum. If the spatial correlation distance decreases, it 
is harder to achieve localization. 

In Fig. [5] wc display dependence of the density 
PB{x,x,t) versus a; for i^ = 8.67 x 10~^ {uj = 10 Hz), 
and three values of (cr,Vb): (0.13,0.465), (0.25,0.478), 
and (0.40, 0.485). Note that Vq can be regarded as a con- 
stant close to 0.47 and we will omit to explicitly write the 
values of Vq besides a in further text; the variations of 
Vo are a consequence of the method utilized to construct 
the random potential. We observe that the exponential 
tails decay faster for larger values of a. 

In order to underpin our observations we compare our 
numerical results to the predictions of the formalism pre- 
sented in Refs. [2a, E^li which has been used to study 
Anderson localization of a Bose-Einstein condensate. In 
the approach of Ref. [2Q|, the wave packet is consid- 
ered to be a superposition of almost independent plane 
waves (fc-components); the wave packet had a high mo- 
mentum cut-off at the inverse healing length (of the con- 
densate). The assumption that the fc-components are 
almost independent means that we should be able to em- 
ploy this formalism here as well, despite of the fact that 
our wave packet is only partially coherent (that is, it is 
in the Tonks-Girardeau regime, rather than being con- 
densed). However, for our wave packets, the momentum 
distribution does not have a clear cut-off but rather de- 
cays smoothly to zero as fc — ^ 00, and we must adopt 
a somewhat different procedure to determine the decay 
rate of the whole wave packet from the decay rate of the 
fc-components. Within the formalism |26l . |42| . the ex- 
ponential decay-rate of every fc'-componcnt is calculated 
by using perturbation theory, and the decay rate 7(fc) is 
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FIG. 4: (color online) First-order correlations in the ini- 
tial state decay algebraically. Shown are the single-particle 
density matrix |pfl(0, a:,0)| (a), and the degree of first-order 
coherence |^fl(0, a;,0)| (b), at t = for the initial state 
corresponding to uj — 10 Hz. The graphs are plotted for 
three values of a as indicated in the legend. Black dot- 
ted lines depict the fitted curves |p_b(0, a::,0)| ~ |x-|^"'^'* and 
\tiB{0,x,0)\^\x\-°-'\ 



given as a scries 7(fc) = X]n>2 7 ^^ increasing orders 
[4^ . The leading contribution 7'^^ to the decay rate ^(k) 
arises from the Born approximation, wherein |26l . |42| . |43| 



^i2)^k)^^f^ 



cr(l-fccr) 



(12) 



for k < a^^ , and zero otherwise: the next order term is 
given in the Appendix, see also [42| . In order to evaluate 
the decay rate of the expanding wave packet from 7 (A:) 
we adopt the following simple procedure. We assume 
that there is an effective high momentum cut-off khmc, 
which we evaluate as follows: For the initial wave pack- 
ets (determined by the trap frequency), we have numer- 
ically calculated the Lyapunov exponents A determin- 
ing localization. For the potential parameter <t = 0.13 
(Vo = 0.465), and the initial condition corresponding to 
a; = 10 Hz (for these parameters the system is close to 
the Born approximation regime (26[), we choose the ef- 
fective high momentum cut off khmc such that it fits the 
numerically calculated decay rate, that is we extract khmc 
from the equation A = ^{khmc)] this yields khmc = 1.79. 
Then, by using this value, we calculate the Lyapunov 
exponents j{khmc) in dependence of a. In Fig. [3] we il- 
lustrate the functional dependence j{khmc) vs. a; dotted 
blue line depicts the Born approximation 7"', and solid 
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FIG. 5: (color online) Correlations in the steady (An- 
derson localized) state. The single-particle density matrix 
\pB{0,x,t)\ [blue line in (a), indicated with the arrow], and 
the degree of first order coherence |/is(0, a:, i)| [blue line in 
(b), indicated with the arrow], at the time 4 s. The param- 
eters used are cr = 0.13 and lu = 10 Hz. Red lines (in both 
panels) depict the single particle density pB{x,x,t). Insets 
enlarge the region where \x\ is small and where correlations 
decay approximately exponentially. For |a;| in the region of 
the density tails (|a::| > Lt), \fiB{0,x,t)\ reaches a plateau. 
See text for details. 



red line depicts 7^^-* -f- 7''^-'. We see that within the pa- 
rameter regime studied here the trend is well described 
with the perturbation approach. Quantitative deviations 
occur because higher order terms of the perturbation the- 
ory are not negligible and should be taken into account. 
Next we focus on correlations contained within the re- 
duced single-particle density matrix pB{x,y,t). Suppose 
that we are interested in the phase correlations between 
the center (at zero) and the rest of the cloud (at some 
X- value); the quantity ps(0,x,i) will decay to zero with 
the increase of jxj even if the field is perfectly coherent 
simply because the density decays to zero on the tails. 
In order to extract solely correlations from the RSPDM, 
we observe the behavior of the quantity [44| 



tJ-Bix,y,t) 



PB{x,y,t) 



yypB{x,x,t)pB{y,y,t)' 



(13) 



which is the degree of first-order coherence 14J] (in optics 
it is sometimes referred to as the complex coherence fac- 
tor [431). In the context of ultracold gases iiBix,y,t) can 
be interpreted as follows: If two narrow slits were made 
at points x and y of the ID Tonks- Girardeau gas, and 
if the gas was allowed to drop from these slits, expand 
and interfere, pb{x, y, t) expresses the modulation depth 
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FIG. 6: (color online) The single-particle states \tlij{x,t)\'^ 
for j — 5, 9, and 13 in the (Anderson localized) steady state. 
The arrow indicates increase of j. The parameters used are 
CT = 0.13 and a; = 10 Hz. The single-particle states for larger 
j (larger in energy) decay slower with the increase of \x\. See 
text for details. 



FIG. 7: (color online) The absolute value of the real and 
imaginary part of Aij{0,x,t) for i — 1 and j = 13, and five 
different realizations of the disordered potential. The param- 
eters used in the simulation are a = 0.13 and w = 10 Hz. For 



sufficiently large [a;| 
text for details. 



Aij{0,x,t) reaches a constant value. See 



of the interference fringes. In this work we investigate 
correlations between the central point of the wave packet 
and the tails: ^j,b(0, x, t). 

In Fig. |4] we show the averages of the magnitudes of 
the one-body density matrix \pB{0,x,t)\, and the degree 
of first order coherence |/xb(0, x,t)\, at time f = for the 
initial state corresponding to a; = 10 Hz and for three 
values of a. From previous studies of the harmonic po- 
tential ground-state (e.g., see Ref. [1^ for the continu- 
ous Tonks-Girardeau gas and [3^ for hard-core bosons 
on the lattice) it follows that in a fairly broad interval 
of x- values, both \pB{0,x,t = 0)| and |^B(0,a;,i = 0)| 
decay approximately as a power law |a;|~'''° with the ex- 
ponent 7o — 0.5 [l^,[33' despite of the fact that the den- 
sity is not homogeneous; the density dependent factors 
multiplying the power law are also known [l^, [S^ . We 
have observed that the initial correlation functions are 
well fitted to the power law: \pb{0,x,0)\ ~ |a;|~°-^* and 
|^B(0,a;,0)| - \x\-°-^^ for a; = 10 Hz (for a; = 5 Hz, 



obtain \pb{0,x,0)\ 



-0.60 



-0.55 



and \fiB(0,x,0)\ 



). The power-law decay of correlations indicates 
presence of quasi long-range order. Apparently, the prop- 
erties of the small random potential do not significantly 
affect the correlations of the initial state for the trap 
strengths w, and disorder parameters used in our sim- 
ulations. This happens because the initial single parti- 
cle states are localized by the trapping potential, rather 
than by disorder (their decay is Gaussian). The effect 
of disorder on these states becomes more significant for 
weaker traps, because the disordered potential becomes 
nonegligiblc in comparison to the harmonic term ly^x^ in 
a broader region of space. In fact, we expect that if one 
keeps the number of particles constant, for sufficiently 
shallow traps, disorder would qualitatively change the 
behavior of the correlations in the initial state, in a simi- 
lar fashion as when the trap is absent. However, probing 
Anderson localization by using transport (i.e., expansion 
of an initially localized wave packet), is perhaps more 



meaningful for tighter initial traps, where the initial wave 
packets are localized by the trap rather than by disorder. 

For very small values of |a;|, and for very large values 
(at the very tails of the wave packet) there are deviations 
from the power law behavior [l^, l38| . The behavior of 
\^^b{0, X, 0)1 at the tails, where \pb{0,x, 0)| starts to grow 
up to some constant value is attributed to the fact that 
higher single-particle states "0^(2;, 0) decay at a slower 
rate with the increase of |a;|, and therefore spatial coher- 
ence increases in the tails (see also the discussion below) . 

After the Tonks-Girardeau gas expands in the disor- 
dered potential and reaches a steady-state, the behav- 
ior of pBiO,x,t) and fj,BiO,x,t) significantly differs from 
that at t = 0. This is shown in Fig. [SJ where we dis- 
play the magnitude of the two functions for a = 0.13 and 
a; = 10 Hz. We observe that \pB{0,x,t)\ exhibits a fairly 
fast exponential decay for small values of |a;|, that is, in 
the region where the density is relatively large [see the 
inset in Fig. [5ja)]. This fast decay slows down up to 
sufficiently large values of x, i.e., |a;| > if, where we ob- 
serve slower exponential decay of |pb(0, x, t)\, which cor- 
responds to the exponentially decaying tails in the single- 
particle density of the localized steady state. Regarding 
the degree of first-order coherence \p,B{0,x,t)\, we find 
that for sufficiently small \x\, it decays exponentially [see 
the inset in Fig. [5{b)]; however, as x approaches the 
region of exponentially decaying tails |a;| > Lt, the expo- 
nential decay of \iiB(fi,x,t)\ slows down until it reaches 
roughly a constant value in the region |a;| > Lt- This 
plateau occurs because single-particle states ipj decay 
slower for larger j values (they are higher in energy and 
momentum), and due to the fact that for sufficiently large 
|a;| , the matrix elements Aij (O, x, t) , which are important 
ingredients in expression ^ for \pB{0,x,t)\, also reach a 
constant value. This is depicted in Figs. |6]and[7l which 
display \^pj{x,t)\'^ for j = 5,9, and 13, and ^143(0,2;,^) 
(real and imaginary part) for five different realizations of 
the disordered potential. We clearly see that Aij{0,x,t) 
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FIG. 8: (color online) Correlations in the steady (Anderson 
localized) state in dependence of the disorder parameter a. 
Shown are the correlations at t = 1450 (= 4 s) for a = 0.13, 
0.25, and 0.40. The arrows indicate increase of a, whereas 
the vertical bars indicate the uncertainty in correlations at 
the plateau values. For a given set of parameters a and u, 
90% of the simulations (one simulation is made for a single 
realization of Vd) fall within the vertical bars. The initial 
trap parameter is a; = 10 Hz. See text for details. 



reaches a constant value (generally complex off the diag- 
onal) , which differs from one realization of the disorder 
to the next; this is connected to the fact that the integral 
/„ dx'ip*{x' ,t)ipj{x' ,t) converges to a constant value for 
sufficiently large x, which is a consequence of the expo- 
nential localization. The fluctuations in Aij(0,x,t) are 
reflected onto the fluctuations of the plateau value of 
|^b(0, x,i)|. We have compared the averages of the ma- 
trix elements Aij (0, x, t) for large x (at the plateau) for all 
values of i and j. They are all within one order of mag- 
nitude with Ai3_i3(0,a;,i) {N = 13) being the largest, 
more specifically, the averages of some of the absolute 
value in our simulations are 1^13,13(0, x,t)\ ~ 0.25 x 10~^, 
\ATj{Q,x,t)\ =0.16x 10-3, |Ai',i(0,a;,i)| = 0.06 x lO^^, 
and 1^1,13(0, a;,i) I = 0.03 x 10"^. Thus, the values of the 
matrix elements to some extent enhance the contribution 
of the highest single-particle states in the correlations 
|/x(0,x,i)|. It is worthy to mention that identical effect 
is observed in incoherent light solitons (e.g., see |39|). 
where the coherence also increases in the tails, which is 
observed in the complex coherence factor in optics (in 
the case of solitons, it is nonlinearity, rather than disor- 
der which keeps the wave packet localized). 

Figure [5] displays the correlations for three values of 
(7 (cr = 0.13, 0.25, and 0.40). In the parameter regime 
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FIG. 9: (color online) Correlations in the steady state in 
dependence of the initial trap parameter u. Shown are the 
correlations at i = 1450 (4 s) for two initial conditions corre- 
sponding to cj = 5 Hz and 10 Hz. The arrows indicate the in- 
crease of (jj, whereas the vertical bars indicate the uncertainty 
in correlations at the plateau values. Other parameters are 
identical as in Fig. [T] 



that we investigated, we found no clear dependence of 
|p_B(0,a;,t)| and |/i_B(0,a;, t)| (in the steady state at 4 s) on 
(T for small values of |a;|. For |a;| in the region of the tails, 
\x\ > Lt, the correlations \^B{^,x,t)\ asymptote larger 
values for larger a. All of the qualitative observations 
above were made throughout the parameter regime that 
we investigated numerically. 

In Figure [9] we display \pB{^,x,t)\ and |^B(0,a;, t)|, for 
two different initial conditions corresponding to oj = 5 Hz 
and 10 Hz; other parameters are identical as in Fig. [1] 
that is, a ~ 0.13. We have observed that for asymptotic 
values of f, the magnitude of correlations |^B(0,x,t)| is 
lower for larger values of w. It is worthy to mention that 
a plateau in |/iB(0,x,t)| occurs in every numerical simu- 
lation for a given realization of the disordered potential, 
but in every realization on a somewhat different value; 
the vertical bars in Figs. [8] and [9] indicate a spread in 
plateau values in our simulations. 

If we compare Fig. [5] with Fig. [5J and Fig. [S] with 
Fig. [1] (compare the direction of arrows in corresponding 
figures), one finds that if the density in the tails decays 
faster, the correlations between the center of the cloud 
and the tails decrease slower. This is in agreement with 
our interpretation that a slower decay (with \x\) of higher 
single-particle states ipj (see Fig. [6]) is responsible for 
the creation of plateaus in |/i(0,a;,i)| and increase of co- 
herence in the tails; namely, if the density decays faster 



(with |a:|). the highest single particle state tpN will be- 
come dominant for smaller values of x leading to greater 
coherence in \^{0,x,t)\. 

Let us now extrapolate our numerical calculations and 
results to larger particle numbers. Suppose that we keep 
all parameters fixed, and increase only N. The energy 
of the initial state as well as the high momentum cut-off 
khcm increase with the increase of N. Our simulations up 
to a finite time up of 4 s would not be able to see exponen- 
tially decaying tails of the asymptotic steady state. By 
employing the results of Ref . ^^ , one concludes that the 
steady state will always be localized, however, at larger 
values of N, the Born-approximation mobility edge [42| 
will be crossed and the exponents describing the expo- 
nentially decaying tails will be smaller. The plateaus in 
the correlations will still exist in the regions of these tails, 
however, the value |/iB(0, x, t)\ will decrease with the in- 
crease of N (simply because more single particle states 
tjjj are needed to describe the Tonks-Girardeau state), 
and both the exponentially decaying tails together with 
the plateaus will be harder to observe. The effect where 
the coherence of the localized steady state increases in 
the tails should however be observable also with partially 
condensed BECs, below the Tonks-Girardeau regime. 



disordered potentials, Anderson localization within the 
framework of the Lieb-Liniger model describing a ID 
Bose gas with finite strength interactions (which becomes 
identical to the Tonks-Girardeau model when the inter- 
action strength becomes infinite). These studies should 
provide further insight into the influence of wave coher- 
ence (within the context of optics), and the influence of 
interactions on Anderson localization (within the context 
of effectively ID ultracold atomic gases). 
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Appendix A: Construction of the disordered 
potential 



IV. CONCLUSION 

We have investigated Anderson localization of a Tonks- 
Girardeau gas in continuous potentials [Vd (x)] with con- 
trolled disorder, by investigating expansion of the gas in 
such potentials; for the initial state we have chosen the 
Tonks-Girardeau ground state in a harmonic trap (with 
Vd (x) superimposed upon it) , and we have analyzed the 
properties of the (asymptotic) steady state obtained dy- 
namically. We have studied the dependence of the Lya- 
punov exponents and correlations on the initial trap pa- 
rameter Lu [5—10 Hz], and the correlation length of the 
disorder a [0.13 — 0.40 /im]. We found that the Lyapunov 
exponents of the steady state, decrease with the increase 
of V. In the parameter regime considered the Lyapunov 
exponents increased with the increase of a, which was un- 
derpinned by the perturbation theory. The behavior of 
the correlations contained in the one-body density ma- 
trix pB{x,y,t) and the degree of first order coherence 
indicate that the off diagonal correlations |p_b(0, x, t)\ de- 
crease exponentially with the increase of |a;|, due to the 
exponential decay of the density, however, in the region of 
the exponentially decaying tails, the degree of first-order 
coherence |//B(0,x,i)| reaches a plateau. This is con- 
nected to the behavior of the single-particle states used 
to construct the Tonks-Girardeau wave function and to 
the increase of coherence in the exponentially decaying 
tails. This effect is analogous to the one found in inco- 
herent optical solitons, for which coherence also increases 
in the tails. 

As a possible direction for further research we envision 
a study of Anderson localization for incoherent light in 



In this section we describe the numerical procedure uti- 
lized for construction of the disordered potential Vd{x). 
The x-space is numerically simulated by using 33000 
equidistant points in the interval x g [—2000,2000]. 
From this array, we have constructed a random array 
V = exp[27rrand(a;)i] of the same length, where rand(.T) 
denotes a random number in between and 1. Then 
we calculated a discrete Fourier transform of v (call it 
v{k)), and introduced a cut-off wave vector Kcut- Vd{x) 
was chosen to be an absolute value of the inverse dis- 
crete Fourier transform of v{k)Q{2k / Kent) [where 6(a::) 
is one for |a;| < 1 and zero otherwise]. We have calcu- 
lated the autocorrelation function Ac [x) of the potential 
and fitted it to the functional form snv{x / a) / {x / a)"^ to 
get the correlation length a. The autocorrelation func- 
tion Ac{x) of the disordered potential Vd{x) is identical 
to the autocorrelation function of the potential used in 
the experiment of Ref. M] , and the theoretical studies 
conducted in Refs. |2a,|42|- The higher order correlators 
differ, but they do not qualitatively change any conclu- 
sions in the parameter regime studied here. 

Let us define the term ^'^'^\k) which is used in the 
series 7(fc) = E„>2 T^"-* (see Ref. [13 and Fig. ^. This 
term depends on the 3-point correlator of the random 
potential ^■. 

C-i[ , ) = —3 {Vd{,x' - Xi)Vd{,x' - X2)Vd{x))x', 



a a- Vi 



and it is given by: 



T^^Hfc) = ^/3(M, 



(Al) 
(A2) 



where /3(k) is defined as: 



dv C3{u,v) sin{2Kv) . (A3) 



The 3-point correlator was calculated numerically to ob- 
tain Fig. El 
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